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METHOD AND APPARATUS FOR REDUCING THE CREST FACTOR 
OF A SIGNAL 

5 

j Background 

The present invention relates to a method and to an apparatus for 
reducing the crest factor (the ratio of peak to average value) of a signal, for 
example, a multi-carrier communication signal. 
10 In recent years multi-carrier communication systems have been widely 

used in particular for x DSL communication systems (Digital Subscriber Line) 
like ADSL (Asymmetric Digital Subscriber Line) or VDSL (Very High Speed 
| Digital Subscriber Line). F - tfc-7-frhows Figure 1 illustrates a schematic block 
diagram of such a transmission system. A serial data signal a is fed to a 
15 | serial/parallel converter 1 whicfa that converts the serial digital data signal a into 
data packets with N/2 sub-packets, N being an even number. One data packet is 
| transmitted in parallel to an encoder 2, which assigns each sub packet to a 

separate carrier frequency and supplies a first digital signal vector to an inverse 
| Fourier transformer 3^ which performs an inverse Fourier transformation on this 
20 vector and generates a second digital signal vector comprising N samples of a 
signal to be sent. This second digital signal vector is transmitted to a 
parallel/serial converter 23, which supplies the elements or samples of the 
second digital signal vector to a digital filter 24 followed by a digital-to-analog 
| converter 25 and a line driver 26. The thus-generated analog transmit signal is 
25 transmitted via a channel 27, whereby noise b is added, symbolized by an adder 
28. On the receiver side, the signal is equalized by an equalizer/an analog-to- 
digital converter 29. Then the signal is decoded by performing the reverse 
operations of the encoding elements 1 to 23, namely through a serial/parallel 
converter 30, a Fourier transformer 3 1 , a decoder 32, a slicer 33 and a 
30 parallel/serial converter 34. 

Such a communication system is for example disclosed in US 6,529,925 
Bl, the contact of which is incorporated by reference herein. 

Since the transmit signal is composed of a plurality of different signals 
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having different carrier frequencies and amplitudes and phases being determined 
by the data signal and thus having no predetermined relationships, the amplitude 
of the transmit signal has approximately a Gaussian distribution. F - ig - . -- 8 
shows Figure 2 illustrates the probability h of the amplitude A of the transmit 
signal as determined by a simulation for a discrete multitone modulated transmit 
signal with a Fourier transform block length of 256. 

Because of this Gaussian distribution, the crest factor of the signal is 
rather large, that is, the transmit signal has a rather high maximum amplitude 
value compared to the effective or average amplitude value. Since both the 
digital-to-analog and analog-to-digital converters^ as well as the line drivers,, 
have to be adapted to handle the whole possible amplitude range, these elements 
have to be defined accordingly causing additional costs and chip space. It is 
therefore desirable to reduce the crest factor, that is to reduce the maximum 
amplitude. 

In principle, two different approaches are known to reduce the crest 

factory 

1t-N 4A first m ethod for reducing the maximum amplitudes whkh-disturbs 

the transmit signal..- These methods comprise clipping methods as described for 
example in US patent no. 6,038,261. 

3r- MA second m ethod for reducing the maximum amplitude without 

disturbing the signal. 

In general, these methods use one or more of the carrier frequencies to 
modify the transmit signal in order to reduce the maximum amplitude. The 
carrier frequencies used for this purpose may »et-or may nolonly partially be 
used for the actual data transmission. 

One of these methods is described in the already cited US patent 
6,529,925 Bl. There, the Nyquist frequency is used as a single carrier frequency 
for correction purposes, that is^ the last frequency in the inverse Fourier 
transform. In an ADSL signal,, this frequency is not used for data transmission so 
that the correction does not influence the transmission capacity. However, the 
performance of this method is limited since only a single carrier frequency is 
used for correction. Furthermore, this method is not applicable to VDSL signals 
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since the Nyquist frequency is outside the usable frequency range both for 
downstream and for upstream transmission. 

In US 6,424,681 Bl a method for reducing the crest factor using a 
plurality of carrier frequencies is disclosed. These carrier frequencies are 
preferably evenly distributed over the whole usable frequency range. From these 
carrier frequencies a normalized correction signal, a so-called kernel, is 
generated which has a "Dirac"-like shape, that is, * tieh % that comprises a 
single peak as far as possible. To correct a transmit signal, this correction signal 
is phase shifted to the position of the maximum of the transmit signal and then 
scaled with a suitable scaling factor depending on the maximum amplitude of the 
transmit signal. Then x this correction signal is subtracted from the transmit 
signal. This can be repeated several times to iteratively correct several maximum 
or peak values. For transmission systems with a great number of carrier 
frequencies and consequently a great number of signal values in each frame, like 
a VDSL transmission system, this method is difficult to realize since it needs a 
relatively long computation time. Furthermore, through the use of a kernel, the 
carrier frequencies used for the correction have to comprise both low and high 
frequencies which, consequently, are not usable for data transmission. The use of 
low carrier frequencies, on the other hand, leads to a greater loss of transmission 
capacity since lower carrier frequencies can be modulated with more bits than 
high carrier frequencies due to the lower damping. 

Summary 

It is therefore an object One embodiment of the present invention te 
provides a method and an apparatus whic ht hat effectively reduce the crest factor 
using a limited number of carrier frequencies. Furthermore, i- t - i - s - an - ebj e e - t - t - e 
pr - ev - ide one embodiment provides such a method and such an apparatus 
which that are usable for VDSL transmission. 

According to one emb< i . \ the invention^ for reducing the crest 

factor of a signal using a plurality of partial correction signals having 
predetermined frequencies, the following steps are carried out: 
(a) determining a time position of a maximum absolute amplitude of the 

signal, 



3 



PATENT 
1435.106.101 
13968US 

(b) calculating an amplitude and a phase for the respective partial correction 
signal depending on said maximum absolute amplitude and said time 
position determined in step (a), 

(c) subtracting the respective partial correction signal from said signal to 

5 obtain a partially corrected signal which is used as the signal in step (a) 

for the next one of the plurality of partial correction signals, and returning 
to step (a) for calculating an amplitude and a phase for the next partial 
correction signal, 

said method further comprising the step of 
10 (d) outputting the last obtained partially corrected signal as the corrected 
signal having the reduced crest factor. 

As for each of the predetermined frequencies, i.e. carrier frequencies, an 
| amplitude and a phase is calculated, and it is possible to use the predetermined 
frequencies available in an optimum manner to correct the signal. 
15 Steps (a) to (c) may be repeated a given number of iterations to obtain 

even better results. 

In step (b), the amplitude may be calculated according to 

A = g ■ (max{X(0 • cos(2^T (t -t max))}+min{x<7) • cos(2;zf (t - 1 max))}) 

20 

A being the amplitude, g being a predetermined factor, / being the respective 
predetermined frequency, t being the time, / max being said time position, X(t) 
being said signal and max and min designate the maximum and minimum 
operator, respectively. The phase accordingly amounts to Inf ■ t max . 

25 For discrete multitone modulation signals as mentioned in the 

introductory portion, the signal may be represented as a signal vector of signal 
values at N sampling times. Accordingly, the above formula may be 
reformulated basically by replacing the time by the number of the sample and 
replacing the frequency by a number of the frequency divided by N . 

30 The method as described so far is suitable for signals where the signal 

vector does not comprise too many samples. For transmission systems like 
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VDSL systems carrying out the described method in full would cost 
considerable calculation time. 

Therefore, instead of performing the above method on the signal or on a 
vector representing the signal, it is preferred to one embodiment performs the 
method on a vector containing only a predetermined number of maximum 
amplitude values of the signal. This predetermined number may be significantly 
lower than the number of samples in the actual signal vector, therefore saving 
considerable calculation time while only marginally lowering the performance. 
To do this, the positions of the elements of the vector with the maximum 
amplitude values in the original signal vector have to be stored since the final 
correction has to be performed on the signal itself. 

Br ief Description of the Drawings 
The accompanying drawings are included to provide a further 
s ' Hiding of the present invention and arc incorporated in and constitute m 

part of this specification , , ' 1 \ illustrate the enu Ui i> re cm 

.mention a n I m i iplion serve to explain ho \ iw-ph • \ 



i « ' ' ^ > ~ i, < , j Me! m m> h 



ntended advantages of the r ji irj \ on will be readily appreciated as they 



become better understood by reference to the following detailed description. 
The cL i i i i > i h tls i 

! isf.ieii'.e iiiiiii>.i j! ■J>>i^i C' ■!! .-iH'iiJii;^ - in ill u p n'\ 

Fi-g.--4--i sFigure 3 illustrates an embodiment of an apparatus for reducing a 
crest factor of a signal according to the present invention-. 

F igs. ? \ 'gas - I A and 54B 4tow iihistratcs simulations of the 
performance of the method of the present invention for ADSL signals^ 

- F - igs - .-3 Figures 5 A and 3-5 B show illustrate further simulations of the 
performance of the method of the present invention for ADSL signals-,-. 

Figs -11 gures 6A and 4 B %h 1 Hi i usable frequency ranges for 
VDSL in a downstream and in an upstream direction, respectively^ 

F%sh- 4Fi gures 7 A and 3-TB illustrate simulation results for VDSL 
upstream^. 
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F i gs - . -- 6 Fi g ures 8 A and 68B show simulation results for VDSL 
downstream-,.. 

Fig.-7-Bhewt J igurc 1 illustrates a standard multi-carrier transmission 
system T ,-an4 

Fife-8-shewsFigure 2 illustrates an amplitude probability distribution for 
a standard multi-carrier transmission system. 

As already described in the introductory portion with reference to Fi-gv 
^Figure 1 , a transmit signal in a multitone transmission like discrete multitone 
transmission comprises a number of samples derived from parallel processing of 
a number of bits of serial data, a data block. This transmit signal may be 
described as a vector 

X T =[x(l),x(2),...,x(AT)L (1) 

N being the number of samples and x(n) being the respective samples, n 
ranging from 1 to N . The index n thus denotes the time position of the 
respective sample. ".7"' indicates that the vector in equation (1) is written in a 
line instead of in a column. 

TOne embodiment of the present invention is to determine a correction 
vector Xk so that the maximum absolute value or amplitude of the elements of 
the vector Xs with 



(2) 



assumes a minimum value. The correction vector Xk is a superposition of 
several partial correction vectors corresponding to a number of carrier 
frequencies or carrier tones reserved for forming the correction vector Xs , i.e. 

Xk = Y^Xk i , (3) 
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Nt being the number of carrier frequencies reserved for correction. Xk i denotes 
the i -th correction vector contribution of carrier frequency number i . 
In general, the components of Xk t may be written as 



xk t (k) = a,Qi)co^27F/i- j + b t (ju) sin^2^u j , (4) 



wherein k is the component or sample index of the vector, ranging from 1 to N , 
and /li is the number of the respective carrier frequency used to form the 
correction vector xk t assuming that all the carrier frequencies used including 

10 those used for transmitting the actual information are numbered consecutively 
and are spaced evenly apart from each other, starting with 0. Such a numbering 
of carrier frequencies is usually used for example for the carrier frequencies of 
ADSL or VDSL transmission. Generally, equation (4) is an oscillation with the 
frequency determined by |a and an amplitude and phase determined by a t and 

15 b, . 

In the following, an iterative procedure for correcting the vector X and 
for determining the correction vector Xk is given, comprising the following 
steps: 

1 . determining the element of the vector X having the maximum 
20 | absolute amplitude \X(k max)| and its position within the vector X k max , and 

2. forming an auxiliary vector Xh according to 



= x(k)cos(2^ k kmaX \k = 1,2,..JV, (5) 



xh being the components of Xh . Since the cosine term for k = k max is equal 
to 1 , the auxiliary vector has the element with the same maximum absolute 
amplitude at the same position k max as the vector X . 

3. For the carrier frequency ju , a partial correction is carried out 
according to: 
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*(k) new = x(k) M - g(max{xh(k)} + mm{xh(k)}) ■ 0.5 • cosj^ - j,.„,„„(6) 
k = \,2,-,N> 

5 wherein the indices "new" and "old" indicate that the elements of the vector X 
are replaced by the new elements. In equation (6), max is the maximum operator 
yielding the maximum value of all the xh(k) and min is the corresponding 
minimum operator. In general, the minimum will be a negative value. It should 
be noted in this respect that the maximum absolute amplitude determined in step 
10 1 may be either the maximum or the minimum value. The factor g is an 

appropriate converging factor which may be chosen to be 1 or may vary from 
iteration to iteration as explained below. The factors 0.5 and g may be drawn 
into a single factor. 

j 4. Ralso given is r epeating steps 1 to 3 for all carrier frequencies ju 

15 used for correcting the signal, whereby the "new" vector X is used for the 

respective next carrier frequency, 
j 5. R also given is r epeating steps 1 to 4 L times. The converging 

factor g can be chosen to decrease from iteration to iteration ensuring a better 

convergence. 

20 The total correction vector Xk would then be a sum of all the corrections 

carried out in step 3. 

The vector X T = [x(l),x(2),...,x(N)] resulting from this method has the 

smallest maximum absolute amplitude possible with a correction signal 

consisting of the given correction carrier frequencies. 
25 On the other hand, for vectors X having a large number of elements this 

algorithm needs a large realization effort since for each iteration for each carrier 

frequency used for correction a correction term has to be subtracted 

componentwise from the original vector X . For example, for a VDSL 

transmission vector the vector X has 8192 elements. 
30 Therefore, a simplification of the above algorithm is needed for vectors 

X having a large number of elements. The general idea for simplifying the 
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above algorithm is not to perform the algorithm on the complete vector X , but 
on an auxiliary vector Xm T = [xm(l),xm(2),...,xm(M)] comprising the M 
elements having the largest absolute amplitudes of the vector X , M being 
much smaller than N . For example, for VDSL systems M may be chosen to be 
5 32 which is considerably smaller than 8192, thus saving considerably 

computation time. Since the correction itself has to be performed on the whole 
signal, that is on the vector X , a further auxiliary vector 
Pm T = [pm(\), pm(2),..., pm(M)] is needed for storing the positions of the 
elements of the vector Xm within the vector X , i.e. xm(k) = x(pm(k)) . 
10 An algorithm for determining the vectors Xm and Pm will be given 

later. 

In the following, it will be shown how the algorithm described above has 
to be performed using the vector Xm . The following steps have to be carried out 
corresponding to the respective steps of the algorithm already explained: 
15 1. the position k max of the element of the vector Xm having the 

| largest absolute amplitude or value \xm{k max)| is determined and 

2. an auxiliary vector Xmh according to 

xmhik) = xm(k) ■ cos(^ pm{k) - ^ max) j; k = 1,2 A/ (7) 

20 

is calculated, wherein the xmh are the components of the vector Xmh . Thus, 
through the use of the vector Pm within the cosine term, the cosine term 
| assumes the "correct" values for the elements of the vector Xm whicht hat 
correspond to the values of the correction finally performed on the whole signal 
25 X. 

3. A partial correction corresponding to the one for the full vector 
X is carried out: 

xm(k) new = xm(k) old - g • (max {xmh(k)} + m\n{xmh(k)}) ■ 0.5 • (8) 
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• cos 2n/j,- 



pm(k) - pm(k max) 



U = 1,2,...,M, 



After performing the algorithm the correction vector Xk for the signal vector, 
X has to be calculated. To this end, it is helpful to store for each partial 
correction and for each frequency ju the correction amplitude 



Wherein / is again the number of the earner frequency as in equation (4) 
and j is the number of iteration. 

4. Repeatin g Finajjv steps 1 to 3 are„ repeated for all carrier 
frequencies ju used for the correction signal or vecto r, and 

5. Repeating steps are repeated 1 to 4 L times, possibly with 
decreasing converging parameter g . 

At the end of this procedure, an auxiliary vector Xm having a minimal 
maximum absolute amplitude is obtained. From the stored amplitude and phase 
values Au(i,j) and Ap(/,y) , the amplitudes and phases for the final correction 
vector Xk with N components can be calculated. Each partial correction vector 
for a single correction carrier frequency ju can be calculated according to 
equation (4). The respective correction amplitudes a ; 0) an( l ^G") can be 
calculated as follows: 



Au(i,j) = g ■ (max{xmh(k)} + mm{xmh(k)}) ■ 0.5 4 



(9) 



and the corresponding phase 



Ap(j,j) = pm(.k im )* 



(10) 




(11) 
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Alternatively, the amplitudes and b t (n) may be computed 

iteratively in step 3 of the algorithm given above, so that the Aw and Ap do not 
have to be stored. In this case, in step 3, the following calculations have to be 
performed: 



The correction vector is composed of cosine and sine values weighed 
with respective amplitude values. The cosine and sine values can be read from a 
sine table or a cosine table. One table is sufficient for the cosine and the sine 
values since its two functions only differ in the phase, that-is the respective 
address of the table read out has to be adapted. The use of such a sine table 
makes the algorithm faster compared to explicitly calculating the sine or cosine 
values each time. 

A reduction of the values to be stored in such a table can be obtained if 
the sine or cosine value is calculated as an interpolation, for example a linear 
interpolation, between stored sine values. It has turned out that the storage of 32 
sine values of a quarter period is sufficient. The remaining three quarters of the 
period of the sine and cosine functions can be calculated using the symmetry of 
these functions. 

It is possible to write the partial correction vector Xk t also as 




(12) 





(13) 



with 



c, (ju) = ^a t 2 (ji) + b t 2 (p) ; <p i (ju) = arctan 



Ami 
U-C")J 



(14) 



11 



PATENT 
1435.106.101 
13968US 



arctan being the arcus tangens operator. 

In this case, only a single sine value has to be calculated or read out from 
the sine table for each partial correction vector with the carrier frequency /j . 
The values for c,.(//) and q> t (ji) from equation (14) may also be calculated 
5 using the known Cordic algorithm. This algorithm is used for calculating 

amplitude of a phase of a complex number when its real end imaginary part is 
given. As a real part a, and as imaginary part b i can be taken. The Cordic 
| algorithm is an iterative algorithm wh i eh that uses only additions and 
subtractions as well as the sign function for determining the sine of a number. 
10 For performing the algorithm, L arcus tangens values have to be stored, L 
being the number of iterations of the Cordic algorithm. After carrying out the 
Cordic algorithm the amplitude of the respective complex number which results 
from the algorithm is enlarged by a fixed factor dependent on L . Therefore, it is 
necessary to divide this value by this factor. To be able to omit these divisions, 
15 the values of the sine table may be divided by the factor in advance. 

A further reduction of computing effort can be obtained if the stored 
values of the sine table, for example 32 values, are multiplied with the respective 
amplitude of the correction carrier frequency in advance and stored in an 
intermediate storage. For computing the partial correction vector this 
20 intermediate storage would simply have to be addressed. No further 
multiplication would be necessary. 

The realization effort can be further reduced significantly by choosing 
the carrier frequencies used for correction in an appropriate manner. If, for the 
representation of equation (13), the carrier frequency is chosen as 

25 

M = 2 t -v (15) 
and the number of elements in the vector X is 
30 N = 2 n (16) 
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which is generally the case for system using inverse fast Fourier transform like 
the system described in the introductory portion, equation (14) transforms to 



xk t (k) = c t (ju) cos^2*2' • v • ^ + <p t (//) j ( 1 7) 



xk, (k) = c, (ju) cos^2;r • v • ^J- + <p M j. (18) 



As can be easily seen, the partial correction vector Xk t is periodic with a 

period of 2"~ l . If the different earner frequencies used for correction differ only 
in the value v , the resulting correction vector Xk only has to be calculated for 
the first 2"~ e values. The whole correction vector is then obtained through 
periodic continuation. 

For the simplified algorithm, as stated above, the auxiliary vector Xm is 
needed containing the M values of the vector X having the largest absolute 
amplitude values. A possible algorithm for obtaining the vector Xm and the 
vector Pm comprises the following steps: 

1. The vector Xm is initialized to contain the M last elements of 
the vector X , ive rthat is, 

xm(k) = x(N - M + k); k = \,2,-M, (19) 

2. The vector Pm is initialized accordingly, fce^ that is, 

pm(k) = N - M + k; k = 1,2,..., M. (20) 

3. A counter X is set to 0: X = 0 . 

4. The element of the vector Xm having the smallest absolute value 
is determined: 

x min = min jjxm (k)\ } (21) 
xmin being that minimum value. 

5. The corresponding position A: min is determined: 

A: min =Position of min|xm(&)|}; ibe rthat is. (22) 

|xm(£min)| = min|x»i(^)|}. 
13 
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6. The counter X is incremented: X = X + L 

7. The element of the vector X designated by the counter A is 
compared with xmin; steps 6 and 7 are repeated until 



5 | \x(X)\ > xmin, 

8. When \x(X)\ > xmin is fulfilled, the minimal element of the 
vector Xm is replaced by the element of the vector X designated by the counter 

| X , and the corresponding element of the vector Pm is replaced by X , jbe rthat is, 
xm(k min) = x(X) 

10 | pm(k min) = X . 

9. The procedure is continued at step 6 until X has reached the 
value N-M 

When this procedure is completed, the vector Xm contains the M 
values having the largest absolute amplitudes of the vector X , and the vector 
15 Pm contains the corresponding positions. 

The time needed for performing this algorithm depends on the arbitrary 
values of the starting vector for Xm . The more large values this vector contains 
when starting the procedure, the less often the contents of the vector Xm has to 
be overwritten and the minimum element of the vector Xm be determined. 
20 Consequently, through a pre-sorting of the vector X it is possible to optimize 
this algorithm. 

j Fig : -4--shw-s -Figure 3 illustrates an apparatus suitable for carrying out the 

method of the present invention corresponding to the algorithms described 
above. A data signal a is supplied to a serial-to-parallel converter 1 and 
25 modulated onto a number of carrier frequencies wherein a predetermined 
number of carrier frequencies are not used for transmitting the data, but for 
building the correction signal as described above. On the thus generated signal, 
an inverse Fourier transform is performed in element 3, and the data is supplied 

| to a parallel-to-serial converter 23 A which is generally used for serially outputting 
30 the corresponding signal vector. Up to this point, the apparatus corresponds to 

I the one already described with reference to Fig. 7 in th e-intredueteay-peft-ien-e-f 
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| t his - pat e nt application. F igure 1 above, that is, in converter 23 the vector X is 
stored. The vector X is transmitted to means 4 for determining the maximum 
amplitude xmax of the elements of the vector X . Comparing means 5 compare 
this maximum value x max with a given reference value xref . If x max is 
5 below xref , switches 6 and 7 are opened, meaning that no correction algorithm 
is performed as the maximum value xmax is below xref which represents a 
maximum tolerable value for the amplitudes or values of the vector X . In this 
case, the vector X is output via a subtracter 35 unchanged, since switch 7 
supplying the negative input of subtractor 35 is opened. 

10 If, however, xmax exceeds xref , the switches 6 and 7 are closed. Via 

switch 6, the vector X is supplied to sorting means 8 which, with the help of the 
parameter M already described above, determines the auxiliary vectors Xm and 
Pm which are supplied to calculating means 9. Calculating means 9 performs 
the above-described iterative algorithm on the vector Xm and then computes the 

15 amplitude and phase values c, and (p, using the frequencies ju allocated for 
correction. At most L iterations are performed. If, however, xmax falls below 
xref before the L iterations are performed, the algorithm is terminated and the 
values for c s and <p t are output. In building means 10, the total correction vector 
Xk is built as described above and supplied via switch 7 to the subtractor 35 

20 where it is subtracted from the vector X . 

It should be noted in this respect that the vector representation serves as a 
means for easily representing the signals. However, the whole procedure may as 
| well be viewed as carried out using the signals themselves, irfe -that is, emitting 
corresponding correction signals having the respective frequencies ju . 

25 It is possible to use less frequencies ju for correcting the signals than 

actually allocated. For example, twelve frequencies may be allocated for 
correction, two or three of them may be used. These two or three frequencies 
should be changed from vector X to vector X , i.e. from frame to frame, to 
distribute the power of the correction signal over all correction frequencies. 

30 Therefore, it is necessary to modify the algorithm slightly. If after the L- 

th iteration xmax is still greater than xref , the procedure is repeated with a 
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different choice of correction frequencies. If this is not the case after a certain 
number of tries, the frequencies which yielded the smallest value xmax are 
used. 

In the following, the performance of the method according to the 
invention will be demonstrated using simulation results. 

Example 1: Downstream transmission in an ADSL system 

For the inverse fast Fourier transform in ADSL systems generally 265 
frequency values which are equally spaced from 0 to half the sampling 
frequency are defined. Therefore, a frame or vector X comprises 512 signal 
values, i-.-e-. thal is, N =5 12. The distance between carrier frequencies is 4.3 125 
KHz, resulting in a sampling frequency of 2.208 MHz. For data transmission the 
frequencies numbers 33 to 255 are used (142.3 to 1 100 KHz). Two different sets 
of parameters were simulated. The first simulation was performed using 
frequency numbers 254, 217, 247, 225, 239, 231, 210 and 243 for correction 
puiposcs. M was set to 8, L, the maximum number of iterations, also to 8. xref 
was set to 4. 1 . The power of the signal was normalized to 1 , so that the peak 
value corresponds to the crest factor. 

■FtgS:-2 -Figures 4 A and -24 B <-b ' II mucs the results for these values. 
Fig. 2 Figure 4A shows iliustrates the probability p for the occurrence of various 
crest factors C given as a ratio, Fi-g.---2-B--i-s -Ftgure 4B illustrates the same graph 
with the crest factor c given in decibel. Curve 1 1 shows the theoretical Gaussian 
distribution. Curve 12 shows the result without correction. The reason why curve 
12 deviates from curve 1 1 is the limited simulation time, for a longer simulation 
time also the even higher crest factors would eventually occur. It can be seen that 
crest factors above 5.5 occur with a probability greater than 10~ 8 . 

Curve 13 shows the results using the method according to the present 
invention. It can be seen that the probability for a crest factor of 4. 1 or 12.25 dB 
is 10" 8 corresponding to a reduction of 2.9 dB compared to the non-corrected 
case of curve 12. 

For a second simulation, only five carrier frequencies were used for 
correction, namely numbers 240, 224, 208, 192 and 176. These five carrier 
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frequencies are evenly spaced apart resulting in a periodic correction signal or 
correction vector Xk with a period of 32. M and L were both set to 8 as in the 
first simulation, and xref was set to 4.3. The result as shew niilustrated in figfc 
^Fi gures 5 A and 35B corresponding to ■Fia»:-2-Figure8 4 A and 2-4B of the first 
5 simulation. Cuive 1 1 again is the theoretical Gaussian distribution, curve 14 is 
| the uncorrected curve corresponding to curve 12 of Figs-.--2 -Figures 4A and 24B, 
and curve 15 is the corrected curve using the method of the present invention. 
The deviations between curves 12 in Figs. 2 Figures 4 A and 2-4B and curves 14 
in Figs. 3 Figures 5A and *5B again stem from the statistical nature of the 
10 amplitude distribution and the limited simulation time. In this case, a crest factor 
of 4.4 corresponding to 12.85 dB occurs with a probability of 10" 8 . Here, still a 
reduction of 2.3 dB compared to the uncorrected cuive is obtained. 

Consequently, it can be seen that the method of the present invention 
leads to a considerable reduction of the crest factor for ADSL transmissions. 

15 

Example 2: VDSL transmission 

In VDSL systems 4096 frequency values are equally spaced apart from 0 
to half the sampling frequency, resulting in a frame or vector X having 8192 
values or elements. The distance between carrier frequencies is 4.3 125 KHz 
20 corresponding to the ADSL value, resulting in a sampling frequency of 35.328 
MHz. 

For downstream and upstream transmission different frequency ranges 
are defined^, which are frhewn iilustrated in F4g*h- 4Figures 6 A and 46B. Figr 
4Figure6A skews -illustrates the frequencies reserved for downstream 
25 transmission corresponding to frequencies number 257 to 695 and 1 182 to 1634. 
F - ig - ."4 Figure 6 B shew &illustrates the frequencies reserved for upstream 
transmission, • in.;; F frequency numbers 696 to 11X1 and 1635 to 2782. 

First, a simulation for upstream transmission was performed. 

Twelve possible carrier frequencies were allocated for correction 
30 purposes, frequencies number 2688, 2624, 2560, 2496, 2432, 2368, 2304, 2240, 
2176, 2112, 2048, 1984. Three of these frequencies were used for actual 
correction. The allocated carrier frequencies are equally spaced apart and have a 
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distance of 64 (i.e. 64 x 4.3125 KHz), resulting in a correction signal having a 
period of 128 independent of the actual choice of the three carrier frequencies 
used for correction. As parameters for the method of the present invention 
M = 32 , xref = 4.3 and L = 8 were used. A maximum number of twelve 
5 choices of the three carrier frequencies were tried for each correction. 

Figs. 5 Figures 7 A and &TB ^h^- \\u± • ij the simulation results, the 
representation of the results again being similar to those of Figs. 2 Figures 4 and 
-35. Curve 1 1 again represents the theoretical Gaussian distribution, curve 16 the 
signal without correction and curve 17 the signal with correction according to 

10 the method of the present invention. A probability of 10~ 8 corresponds to a crest 
factor of 4.5 or 13 dB according to curve 17, which again is a considerable 
reduction compared to the 5.6 or 15 dB of the uncorrected curve 16. 

For a downstream simulation, six carrier frequencies were used, namely 
frequencies number 1600, 1536, 1472, 1408, 1344 and 1280. The distance 

15 between the carrier frequencies again is 64, which again results in a periodic 
correction signal having a period of 128. For the simulation, the parameters 
| M = 32, xref = 4.3 and L = 8 were used. F t gs : -6 8A and 68B she willustrate the 
results of this simulation. Curve 18 is the result without correction, curve 19 is 
the result with the correction according to the present invention. A probability of 

20 10" 8 corresponds to a crest factor of 4.65 or 13.4 dB, again yielding a significant 
improvement compared to the uncorrected signal. 

It should be noted that the simulation examples given above are only 
given as an illustration, and other parameters may be used depending on the 
amount of crest factor reduction needed and the amount of computation time 

25 possible. For example, a larger value of M generally leads to a better reduction 
of the crest factor, but needs more calculation time. Other carrier frequencies 
than the one used in the simulations may be allocated for correction purposes. 

Although exemplary embodiments of the invention are described above 
in detail, this does not limit the scope of the invention, which can be practiced in 

30 a variety of embodiments. 
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